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I give an introduction to the effective field theory description of top pair 
production at threshold in e + e~ annihilation. The impact of the summa- 
tion of logarithms of the top quark velocity including most recent results 
at next-to-next-to-leading logarithmic order is discussed. 

1. Introduction 

The so-called "threshold scan" of the total cross section lineshape of top 
pair production constitutes a major part of the top quark physics program 
at a future e + e~ collider. In the Standard Model the top quark width 
Tt ~ 1.5 GeV is much larger than the typical hadronization energy Aqcd- 
In contrast to the J ftp or the T region, it is therefore expected that non- 
perturbative effects are strongly suppressed, and that the threshold line- 
shape is a smooth function of the c. m. energy. Throughout this talk I will 
therefore neglect non-perturbative effects associated with the hadronization 
scale Aqcd • From the location of the rise of the cross section a precise mea- 
surement of the top quark mass will be possible, while from the shape and 
the normalization of the cross section one can extract the top quark Yukawa 
coupling y t (for a light Higgs), the top width or the strong coupling. In the 
past, numerous studies have been carried out to assess the feasibility and 
precision for extracting various top quark properties from a threshold run 
(see e.g. Ref. j^j for a recent study). 

In the threshold region, ^/s ~ 2mt ± 10 GeV, the top quarks move with 
nonrelativistic velocity in the cm. frame. Let us define mtv 2 = yfs — 2mt- 
We see that parametrically \v\ < a s . Because in the loop expansion one 
encounters terms proportional to (a s /v) n (n = 0, 1, . . .) in the amplitude, it 
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is necessary to count a s /v of order one and to carry out an expansion in ot s 
as well as in v. Schematically one needs an expansion of o t i = a(e + e~ — ► tt) 
of the form 



R = ^=^E -)xjl(LO); 



a s , v (NLO) ; a?, a s v, v l (NNLO) j . (1) 



The indicated terms are of leading order (LO), next-to-leading order (NLO), 
and next-to-next-to-leading order (NNLO). I call the expansion scheme in 
Eq. Q fixed-order "perturbation theory , although it involves summations of 
the terms proportional to (a s /v) n . It can be implemented systematically us- 
ing the factorization properties of non-relativistic QCD (NRQCD) [Hj. The 
NNLO QCD corrections to the total cross section were calculated some time 
ago m Ref. d 03 El QDj ■ Surprisingly, the corrections were found to 
be as large as the next-to-leading order (NLO) QCD corrections, and from 
the residual scale dependence in the NNLO result, the normalization of the 
cross section was estimated to have at least 20% theoretical uncertainty. 
It was concluded that the top quark mass in a threshold mass scheme can 
be determined with a precision of 200 MeV or better from the shape of the 
cross section. However, the large NNLO QCD corrections to the nor- 
malization of the cross section jeopardized competitive measurements of the 
top width, strong top coupling, or the top Yukawa coupling. Moreover, the 
large NNLO corrections seemed to indicate that, despite the perturbative 
nature of the tt system, high precision computations might not be feasible. 
The determination of higher orders would be necessary to clarify the feasibil- 
ity of the fixed-order approach. Figure^ shows the vector-current-induced 
cross section a(e + e~ — ► 7* — ► tt) at LO, NLO and NNLO in fixed-order 
perturbation theory for typical choices of parameters and renormalization 
scales. (Top threshold production mediated by the axial vector current is 
suppressed by v 2 and small and not relevant here for the discussions on the 
convergence of the perturbative expansion. See e.g. Refs. [TT1 IT2*] .^ 

One way to understand the large scale uncertainties and the large size 
of the fixed-order NNLO QCD corrections is to recall that the perturbative 
tt threshold dynamics is governed by vastly different energy scales, the top 
mass (m; ~ 175 GeV), the top three-momentum (p ~ mv ~ 25GeV) and 
the top kinetic or potential energy (E ~ mv 2 ~ 4 GeV). This hierarchy of 
scales is the basis of NRQCD factorization j3j , which separates hard and 
non-relativistic effects, and which can be implemented in the fixed-order 
scheme shown above. However, in fixed-order perturbation theory NRQCD 
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matrix elements involve logarithmic terms such as 




hi 





(2) 



which cannot be rendered small for a single choice of the renormalization 
scale \i. For example, since E ~ 4GeV, one finds a s {mt) ln(mf /E 2 ) ~ 0.8 
for [i = nit which is of order unity, and fixed-order perturbation theory 
becomes unreliable. Moreover, at higher orders, fixed-order perturbation 
theory cannot distinguish the scale at which to evaluate a s . Mistaking an 
a s {mt) for an a s (rritv 2 ) is a difference of a factor of two. Both problems 
cannot be addressed systematically within the framework of fixed-order per- 
turbation theory. But they can be avoided in a framework that allows for 
renormalization group improved perturbative computations, where all loga- 
rithmic terms are summed through renormalization group equations. Thus 
a better expansion should have the schematic form 



where the indicated terms are of leading logarithmic (LL), next-to- leading 
logarithmic (NLL), and next-to-next-to-leading logarithmic (NNLL) order. 
I have used the expression "better expansion" because in QCD compu- 
tations one should favor expansion schemes where logarithmic terms are 
summed into coefficients and not contained in isolated form in matrix ele- 
ments. This is even true in cases where the non-logarithmic corrections are 
sizeable too, because possible cancellations between large logarithmic and 
large non-logarithmic terms can be unphysical. To accomplish a renormal- 
ization group improved calculation one needs to employ a more sophisticated 
effective field theory approach than the one that is used for the fixed-order 
computations. 

In this talk I give a basic introduction to "velocity NRQCD" (vN- 
RQCD) Ref. Q31 [H EH [TJ] (for a review see also Ref. [Toj), which is an 
effective theory for heavy non-relativistic quark pairs that, through renor- 
malization, sums all logarithms involving ratios of the scales mt, |p| and 
E. It is assumed that E 3> Aqcd- The program is similar in spirit to 
let's say summing QCD logarithms of Mw/rrib for the electroweak Hamil- 
tonian describing b-quark decays. However, in the non-relativistic case the 
program is more complicated because the structure of the relevant degrees 
of freedom in the effective theory action is more involved and because the 
two low energy scales p ~ mv and E ~ mv 2 are correlated through the 
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heavy quark equation of motion E = p 2 /m. At first sight it seems impos- 
sible to keep this correlation of scales and to render the logarithms in the 
effective theory matrix elements shown in Eq. @ small (which is equivalent 
to saying that all logarithms are summed into Wilson coefficients of the 
effective theory). In the effective theory this problem is dealt with by using 
two renormalization group scales in the effective Lagrangian: fig for soft 
(~ mv) fluctuations and for ultrasoft (~ mv 2 ) fluctuations. The matrix 
elements of the effective theory then only contain the logarithms ln(/x|./p 2 ) 
and ln(/z^/ E 2 ). Both renormalization scales are related by flu = fj/g/mt- It 
is therefore convenient to define mv = /xg and mv 2 = fj,jj, where v is a di- 
mensionless parameter. All renormalization group equations of the effective 
theory are expressed in terms of v. Running from v = 1 to v ~ v, where v is 
the typical top quark velocity, sums all logarithms of v and minimizes both 
ln(/i|/p 2 ) and ln(/i 2 y/£' 2 ) in the matrix elements. (Logarithms involving 
rrit are minimized by matching QCD onto the effective theory at v = 1.) 
The correlation of both renormalization scales is essential for the correct 
summation of logarithms. |17j . 

In dimensional regularization the factors of fi e s and multiplying each 
operator in the effective theory action are uniquely determined from its 
mass dimension and v power counting. In this way the scheme indicated in 
Eq. (j3J) can indeed be achieved. 

With this setup and after having formulated the effective theory ac- 
tion the program is similar the one for summing QCD logarithms for the 
electroweak Hamiltonian: 

1. Matching computation of the coefficients Ci(v) of the effective theory 
operators (including external currents) at v = I in a perturbation 
series in a s (mt). 

2. Computation of the anomalous dimensions of the operators of the 
effective action and scaling of Ci{v) from v = \ to v = vq ~ -y~ Cpa s 
using the renormalization group. 

3. Computation of the cross section using the effective Lagrangian and 
the currents renormalized at the low scale v = vq. 

In this talk I show how the vNRQCD effective Lagrangian is constructed, 
and I also spend some time on the points 1.-3. described above including 
new recent results for the total top pair production cross section at threshold 
in e + e~ annihilation. At the end I will briefly discuss the status of top pair 
production at threshold. 
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2. The vNRQCD Lagrangian 

The physical system we wish to describe is that of a heavy quark and 
antiquark with mass m, and energies E ~ mv 2 , and momenta p ~ mv in 
the cm. system where v <C 1. The degrees of freedom from which we have 
to build the effective action can be identified from the relevant momentum 
regions that can be found in an asymptotic expansion of non-relativistic 
scattering diagrams in the cm. frame. ^S] The regions include hard modes 
with momenta (kQ,~k) ~ (m,m), soft modes with momenta ~ (mv,mv), 
potential modes with momenta ~ {mv 2 ,mv) and ultrasoft modes with mo- 
menta ~ (mv 2 ,mv 2 ). Fluctuations with off-shell momenta are integrated 
out because they do not resonate, while for modes that can resonate fields 
are introduced from which we build the effective Lagrangian. Resonating 
heavy quarks can only live in the potential regime, while massless modes 
can live either in the soft or the ultrasoft regime. Thus, the effective La- 
grangian is build from heavy potential quarks and antiquarks (ip p , Xp)> 
soft gluons, ghosts, and massless quarks (A 1 ^, c q , ip q ) and ultrasoft gluons, 
ghosts, and massless quarks (A 11 , c, ip us ). The ultrasoft gluons are the 
gauge partners of momenta ~ mv 2 , while soft gluons are the gauge partners 
of momenta ~ mv. This classification is, of course, only meaningful in the 
cm. frame for a non-relativistic heavy quark-antiquark pair. Double count- 
ing is avoided since ultrasoft gluons reproduce only the physical gluon poles 
where k° ~ k ~ mv 2 , while soft gluons only have poles with k° ~ k ~ mv. 
It is essential that both soft and ultrasoft gluons are included at all scales 
below m because the heavy quark equation of motion correlates the soft and 
ultrasoft scales. We will also see later in Sec.|31 when I show new results 
for the 3-loop (NNLL) anomalous dimension for the ti production current, 
that both soft and ultrasoft running can in general feed into anomalous di- 
mensions induced by potential loops. The dependences on soft energies and 
momenta of potential and soft fields appear as labels on the fields, while 
only the lowest-energy ultrasoft fluctuations are associated by an explicit 
coordinate dependence. Formally this is achieved by a phase redefinition 
for the potential and soft fields 

<Mz)^5>-^ fc (*), (4) 

k 

where k denotes momenta ~ mv and d^(f)k{x) ~ mv 2 (f>k(x). This means that 
the effective Lagrangian contains sums over fields with soft indices, which 
also build up potential and soft loop integrations when we renormalize the 
theory or compute matrix elements. 

The effective vNRQCD Lagrangian for a ti angular momentum S-wave 
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and color singlet state has terms I14[ I15j 
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+E -p^f + ■••- E /4 £ np>p') <^pxl p ^-p 
p p,p' 

+ E 7^%/^ C (P - PO-^A^^T^pxLpT^-p] + • • .0) 
p,p' lP PJ 

where color and spin indices have been suppressed and g s = g s (mtv), g u = 
gsirntv 2 )- All coefficients are functions of the renormalization parameter 
and all explicit soft momentum labels are summed. The covariant derivative 
in the first line contains only the ultrasoft gluon field. There are 4-quark 
potential-like interactions of the form (k = (p — p')) 

v( f\ v c V r (p 2 + P /2 ) v 2 v s 

y(P ' P) = P + 2m? k 2 + rnI + ^f S ' (6) 

where S is the total ti spin operator. Note that the momentum structure 
of the operators satisfies the on-shell conditions and hermiticity. There is 
still some freedom in the choice of the operator basis for the potential in- 
teractions shown in Eq. (jBJ), which can affect the matching conditions and 
the anomalous dimensions. In Refs. ^3 HO] we also used a potential 
interaction term of the form l/(m|k|), while in Ref. [U the order l/(m|k|) 
potentials were implemented only by potential-like interactions with addi- 
tional sums over intermediate indices. 

At NNLL order for the total cross section the coefficient V c of the 1/k 2 
potential has to be matched at two loops [Hj because it contributes at 
the LL level, whereas the coefficients of the order 1/m 2 potentials have to 
matched at the Born level |14j . The l/(m|k|)-type potentials are of order a 2 
and have to be matched at two loops. [2U] There are also 4-quark interactions 
with the radiation of an ultrasoft gluon (last line) and interactions between 
quarks and soft gluons (second line). Due to momentum conservation at 
least two soft gluons are required. The potential terms shown in Eq. @, 
and time-ordered products of soft interactions contribute to the potentials 
that describe the instantaneous interactions between the top quarks. For 
example, the time-ordered product of two soft interactions at leading order 
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leads to the potential interaction [22] 




where the /3j are the coefficients of the QCD beta-function and the a% were 
determined some time ago in Ref . |23| . The sum of Eq. (J7J) and the first term 
in Eq. © agrees with the static potential [21] 1 and constitutes the complete 
1/k 2 potential at NNLL order. The Lagrangian also contains more compli- 
cated interactions describing 4-quark potential-like interactions and 4-quark 
interactions with soft or ultrasoft gluons that contain additional sums over 
intermediate indices. Their form is partly fixed by the renormalizability 
of the theory |2U| . but there is still some freedom in the choice of their 
form. In this talk I use the conventions of Ref. |2J were all l/(m|k|)-type 
potentials are represented by such 4-quark sum operators. I do not display 
their expressions here and refer to Ref. [2^ where they are collected in an 
appendix. 

One might wonder why the effective theory contains both soft and ul- 
trasoft degrees of freedom at the same time, since, naively, one might argue 
that for ti<l the soft fields fluctuate at a much smaller length scale than 
the ultrasoft fields. That both soft and ultrasoft degrees of freedom have 
to be present at the same time is at the heart of the effective theory con- 
struction and is related to the fact that the dispersion relation of the heavy 
quarks, E = p 2 /m introduces a correlation between the soft and ultrasoft 
scales. The presence of soft and ultrasoft degrees of freedom also allows for 
the simultaneous summation of logarithms of p and E, which I mentioned 
in the introduction, after the theory is renormalized. I will come back to 
this point in the next section. 

An important point to mention is that the factors of /ifj and /i% multiply- 
ing the operators in the effective Lagrangian are uniquely determined by the 
mass dimension and the v power counting in d = 4— 2e dimensions |15j . Each 
field in the action is assigned a certain scaling with v to keep its kinetic term 
of order v°. This results for example in ift p ~ (mu) 3//2_<E , ~ (mt) 2 ) 1_f , 
and ~ (mv) 1-6 . Then, since -D M ~ mv 2 , the renormalized combination 
g u -^ must be multiplied by ~ (mv 2 ) 6 for v ~ v so that this gluon term 
also scales consistently as mv 2 . Applying the same procedure to all inter- 
actions results in the factors of fx\j and fi £ s shown for example in Eq. (JSJ). 

1 The statement made in Refs. 1251 [To] that the vNRQCD 1/k 2 potential at NNLL 
order differs from the static potential has shown to be incorrect in Ref. |20|. 
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Note that this procedure is an integral part of the theory and automatically 
leads to the correlation of the two renormalization scales, \x\j oc fj/g/m. 

To incorporate the effect of the large top quark width we include in the 
effective Lagrangian the operators 

^ = E^r t v P + E4^ Xp , (8) 

p p 

where Tt is the total on-shell top quark width. The typical energy of the top 
quarks at threshold is E ~ mv 2 ~ 4 GeV, and in the Standard Model one 
has Tt = 1.43 GeV ~ E. Thus, the propagator for a single top (or antitop) 
with momentum (p°,p) is 

i 

p° -p 2 /(2m) + iT t /2 + ie ' 

which gives a consistent NLO treatment of electroweak effects. 25i The com- 
plete NNLL order treatment of electroweak effects for the total cross section 
is currently unknown. Although it is not expected that the unknown elec- 
troweak corrections are beyond the few percent level this has to be kept 
in mind when I discuss the NNLL QCD corrections. 

Besides the interactions contained in the effective Lagrangian that de- 
scribe the dynamics of the ti pair we also need external currents that de- 
scribe the production of the top quarks. For e + e~ annihilation these cur- 
rents are induced by the exchange of a virtual photon or a Z boson. At 
NNLL order we need the vector S-wave currents J p = ci(v)0pl + c 2 (v)0p2, 
where 

0pl = ^a(ia 2 ) X - p , (10) 
0p2 = — ^p 2 (r(ia 2 )xl p , 
and the axial- vector P-wave current Jp = c%(v) 0p3, where 

0P3 = ^-^[<T,<j-p](ia 2 ) X t p . (11) 

Note that the currents have a soft momentum index. The currents 0p2 and 
0p3 lead to contributions in the total cross section that are t; 2 -suppressed 
with respect to those of the current 0pl. Thus, at NNLL order, two-loop 
matching is needed for c\ and Born level matching for c 2 and C3. Note that 
the two-loop matching result for c\{y = 1) depends on the choice of the 
operator basis used in the effective Lagrangian. 
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3. Anomalous Dimensions and Renormalization Group Scaling 

To sum the logarithmic terms (a s In v) % at NNLL order, as indicated 
schematically in Eq. ©, one needs to determine the anomalous dimensions 
of the potentials and currents that contribute to the cross section at the 
appropriate order. The computation of the anomalous dimensions in the 
effective theory works just as for the effective weak Hamiltonian. The coef- 
ficient of the 1/k 2 (Coulomb) potential, V c , needs to be determined at the 
three-loop level since it contributes already at leading order. The results 
are available in Ref. |2Qj. The anomalous dimensions of the l/rrtf potentials, 
V r ,2,s only need to be determined at one-loop because the 1/mf potentials 
are suppressed by v 2 and they are only needed at LL order. The results have 
been given in Refs. [HI 1217] . (See also Refs. jHUHl-) For the l/(m|k|) po- 
tentials, on the other hand, the anomalous dimensions have to determined 
at two-loops. Also this work has been achieved. [IHJj F° r the current 0pl one 
needs to compute the anomalous dimension at three loops and for 0p2 and 
0p3, they are needed at the one-loop level. The LL anomalous dimensions 
of 0p2 and 0p3 were computed in Refs. |121 119j . I will not give any details 
on these results in this talk since they are already available in the literature 
for some time. 

Here I want to discuss some aspects of the anomalous dimension of the 
current 0pl, which is currently only fully known to NLL order |131 1271 00] . 
but for which I have recently obtained a partial NNLL order result. The 
NNLL order results I will discuss are also of conceptual interest as far as the 
construction and the renormalizability of the effective theory is concerned. 

To obtain the anomalous dimension of 0pl, one needs to determine the 
renormalization constant of the current. Using dimensional regularization 
and the MS scheme one can define the unrenormalized Wilson coefficient as 



At LL order there are no UV divergences that have to be absorbed into the 
current. This means that at LL order the Wilson coefficient c\ does not 
run. 

To determine the N fc LL order renormalization constant of the current 
0pl one has to compute the overall UV divergences of quark-antiquark-to- 
vacuum on-shell matrix elements of spin-triplet S-wave currents at 0(a(? +1 ). 
All lower order UV sub divergences have to be subtracted by lower order 
counterterms. It is mandatory to use matrix elements where the external 



c x = c\ + Set = Z Cl c\ 



(12) 



and one can write the renormalization constant as 




(13) 



10 




(e) 

Fig. 1. Order a 2 s vertex diagrams diagrams for the computation of the NLL anoma- 
lous dimension of c\. The dot and box represent single insertions of the LL 1/k 2 
and the l/mf suppressed potentials, respectively. The circled cross denotes a single 
insertion of the l/(m|k|)-type potentials. Insertions of these potential involve only 
a one-loop diagram because the coefficients of the l/(rn|k|)-type potentials are of 
order a 2 s . 



quarks are on-shell because for off-shell quarks one can obtain UV diver- 
gences which do not belong to the current. This is because we have defined 
our operator basis in the on-shell limit, particularly for the potentials. Most 
of the time one does not have to worry about this issue, but such off-shell 
UV divergences in fact exist at NNLL order. 

Let me first review the NLL result which was first obtained in Ref. |13| . 
The relevant order a 2 s vertex diagrams are shown in Fig.^ There are no 
one-loop sub divergences that have to be subtracted. The computation is 
complicated technically by the fact that there are IR-divergent Coulomb 
phases for on-shell external quarks which have to be distinguished from the 
UV divergences that go into the renormalization of the current. At NLL 
order this distinction can be carried out easily because only 1/e singularities 
occur. From the UV divergences one obtains 



5z, 



NLL 

Cl 647T 2 

2 



+ v^V) + vW(^) + s 2 vj s )(z 



+ a s (wi!/) 



+ aAmv) 



, (14) 



where S 2 = 2 is the squared quark total spin operator for the spin-triplet 
configuration. The terms in the second line arise from a single insertion 
of the l/(m|k|)-type potentials. Using that is renormalization group 
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invariant, one can derive the NLL order anomalous dimension of ci, 



To solve the NLL renormalization group equation one needs to know the 
LL solutions for all the couplings that appear on the RHS of Eq. (|15|) (i.e. 
which mix into c\). From the fact that all higher order l/e n (n = 1, 2, . . .) 
divergences have to cancel in the anomalous dimension (see e.g. Ref. |28j). 
one can also obtain the NNLL order 1/e 2 coefficient 5z^ NLL ' 2 of the renor- 
malization constant Z Cl . [2^ 

The determination of the NNLL order anomalous dimension proceeds 
along the same lines. One can distinguish between contributions from two 
classes. The first class, called "mixing contributions", arises from the terms 
in Eq. (|14|) due to the NLL order running of the couplings on the RHS of 
Eq. (|15|> . The second class involves the computation of three- loop vertex 
diagrams with potential loops and either soft or ultrasoft loops that require 
new ci counterterms. I call this second class "non- mixing contributions". 
It leads to genuinely new contributions in the anomalous dimension of c\. 
By power counting there are no contributions from diagrams with three 
potential loops or which have both soft and ultrasoft loops. I have deter- 
mined the non-mixing contributions recently in Ref. [U and I will discuss 
the computation and the results in some detail in the following. 

At NNLL order the determination of the UV divergences from on-shell 
three-loop vertex diagrams is complicated because there are overlapping 
IR and UV divergences which are difficult to separate from each other. 
An efficient way to avoid these complication is to consider 4-loop current 
correlator graphs rather than the vertex diagrams. The correlator graphs 
are obtained from closing the external quark lines of the vertex diagrams 
with an additional insertion of the current. In Fig. [2] 1 have exemplarily 
displayed the four- loop correlator diagrams with an ultrasoft gluon (loop). 
The imaginary part of the correlator graphs is proportional to the squared 
matrix elements from which all IR-divergent Coulomb phases drop out auto- 
matically. The three-loop (NNLL) renormalization constant of the current 
is then obtained from the three-loop sub divergences of the correlator dia- 
grams that remain after the one- and two-loop sub divergences have been 
subtracted. For dimensional reasons there are in fact no four-loop overall 
divergences in dimension regularization. 

Due to lack of space I cannot give all the details of the computation 
and the full result here and refer to Ref. [2*T]. But I would like to discuss 
some interesting conceptual aspects of the results. First, the NNLL order 
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(a) 



(b) 



(c) 



OIO 



ceo 



(d) 



(e) ,«5«r^ 



OD 



(f) -StWJ^ 




(h) 



(ki 



(i) 




(1) 



(m) 5 s5ffS"^«%. 

Fig. 2. Four- loop graphs with an ultrasoft gluon for the calculation of the ultrasoft 
non-mixing contributions of the NNLL anomalous dimension of c\. The ultrasoft 
gluon is attached to the top quark or to potential-like 4-quark interactions. 



1/e 2 coefficient 5z t 



NNLL, 2 

ci 



I have obtained by the computations agrees with 



the NLL prediction based on renormalization group invariance mentioned 
above. This is a non-trivial check for the consistency of the effective theory 
under renormalization, since it requires that the interactions between the 
relevant degrees of freedom are encoded correctly in the structure of the 
effective Lagrangian. This shows in particular, that the renormalizability 
of the effective theory requires that soft and ultrasoft degrees of freedom 
have to be present for all scales below the quark mass rat. In addition, the 
results for the diagrams in Fig. [21 which involve the soft and the ultrasoft 
renormalization scales, \is and show that the correlation = /x|/mt 
I have mentioned previously is required by the renormalizability of the ef- 
fective theory. As an example, let us consider the contribution to 5Z^ NLL ' 2 
induced by the diagram Fig.|2^ which, when the correlation between us and 
fiu is neglected, has the form 



a s (m t v) 2 a s (m t v 2 ) C\C\ 



TT 



1 



12e' 



1 

6e 



M2 + I-ln(^) 
6 V /4 / 



(16) 



The dependence on In /is and In^u in the 1/e piece that contributes to 
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m = 175 GeV 


m = 4.8 GeV 


V 




U \ v ) 






1.0 


0.0000 


0.0000 


0.0000 


0.0000 


0.8 


0.0069 


0.0041 


0.0308 


0.0425 


0.6 


0.0157 


0.0104 


0.0712 


0.1304 


0.4 


0.0274 


0.0216 


0.1335 


0.4537 


0.2 


0.0435 


0.0512 







Table 1. Numerical values for £ NLL (i/) and £nm LL (^)- The values for m are pole 
masses. The numbers are obtained by evaluation of the analytic results using four- 
loop running for a s and taking ai 5) (175GeV) = 0.107 and ai 4) (4.8GeV) = 0.216 
as input. 



the NNLL anomalous dimension of c\ vanishes only if the correlation = 
(J^/rn is accounted for. This shows that the correlation of the soft and 
ultrasoft scales in the effective theory is required by physical reasons and 
not just imposed by hand. 

Although the computation of the NNLL anomalous dimension of c\ is 
not completed yet, it is instructive to compare the numerical size of the 
non- mixing NNLL contributions with the NLL ones in the running of c\. 
Let us write the solution of Eq. Q15|) for v < 1 as 



In 



ClO) 



ci 1 



C NLL (-) + (^ NLL H + e LL M ) + .... (17) 



where £j^m LL (£m NLL ) refers to the NNLL non-mixing (mixing) contribu- 
tions. In Tab. ^ the values for £ NLL (z^) and f™ L (i/) are displayed for dif- 
ferent v for the top and the bottom quarks. For top quarks we find that 
the NNLL non-mixing contributions are of the same size as the NLL terms 
for the relevant region v ~ v ~ 0.2. Here, the new NNLL order corrections 
shift c\ by about +5%, which is substantial considering that the total cross 
section contains cf. The shift is dominated by the contributions coming 
from the ultrasoft diagrams shown in Fig. [21 However, it is not yet possi- 
ble to draw definite phenomenological conclusions for the normalization of 
the top threshold cross section in e + e~ collisions from this result, because 
the yet unknown NNLL mixing corrections could be sizeable as well. It is 
therefore an important future task to determine the mixing contributions 
to the NNLL anomalous dimension of c\. 

For bottom quarks the NNLL non-mixing contributions are several times 
larger than the NLL terms for the relevant region v ~ v ~ 0.3-0.4. This 
is not unexpected because for bottomonium systems the binding energy 
~ mv 2 is already of order Aqcd- Thus our result seems to affirm that for 
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bb states non-perturbative effects have a rather strong influence, and that 
the vNRQCD description ceases to work even for the ground state. ^2] 



4. Determination of the Total Cross Section 

The total cross section for e + e~ — > 7*,Z* — ► ti at threshold at NNLL 
order has the form 



Otot 0) 



Ana 1 
3s 



F v (s)R v (s)+F a (s)R a (s) 



(18) 



where F v,a are trivial functions of the electric charges and the isospin of the 
electron and the top quark and of the weak mixing angle. At NNLL order 
the vector and axial-vector i?-ratios have the form 

47r r 

R v (s) = — Im c\(v)Ai{v,m,v) +2ci(v)c 2 {v)A 2 {v,m t ,v) , (19) 

S L J 



R a (s) 



47T 



Im 



(20) 



where the time-ordered products of the effective theory currents read (q 
(v^-2mt,0)) 



p,p 



4 x e^ x 







T0pl(x)OL 1 (O) 



, 



A = ~ £ /^e^(o|r[0pl(x)O^ 2 (O) + P 2(x)O^ 1 (O) 



p,p 



p,p' 



0), 
(21) 



The correlators .Aj can be written in a compact form in terms of non- 
relativistic zero-distance Greens functions. For the convention for the l/(m|k|)- 
type potentials employed in Ref. |2J they have the form (v = ((y/s — 2mt + 
iT t )/m t y/ 2 ), 

Ai(v,m,u) = 6N c [G c (v,m,u) + (y 2 (y) + 2V s (i/)) SG s (v,m,u) 
+ V r (v) 5G r (v, m, v) + 5G kin (v, m, v) 



C 2 

C A C F o? s {mv) 5Gc ACF (v,m,u) + -^a 2 s (mv) 5GQ ¥2 (v,m,v 
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+ a 2 (mu)V kl (u) 5G kl (v, m, v) + a 2 s {mv)V k2 {v) 5G k2 (v, m,v) , 
A2(v,m,v) = v 2 Ai(v, m, v) , Az{v,m,v) = — ^- G 1 (a, v , m, v) . (22) 



The Greens function G c contains the sum of 1/k 2 potentials from Eqs. (JJJ) 
and © and was computed numerically in Ref. (See also Ref. 

All other Greens functions were obtained analytically in dimensional reg- 
ularization. The terms 5G S ' r arise from a single insertions of the 1/m 2 
potentials in Eq. © and <5G kin from an insertion of the kinetic energy 
p 4 /(8m|). |12| The terms SG k ' kl,k2 come from a single insertion of the 
l/(m|k|)-type potentials and their expressions were obtained in Refs. |201 
121) . The P-wave Greens function G 1 was also obtained in Ref. ^2] • The 
Greens functions contain UV sub divergences which are removed by the 
renormalization constant of the current 0pl. They still contain overall 
divergences which, however, are not contained in their absorptive part for 
stable quarks. 

I have displayed the result in the pole mass scheme. Since the pole 
mass is plagued by a renormalon ambiguity of order Aqcd , which causes an 
instability in the prediction of the cross section, it is mandatory to switch to 
a "threshold mass" mass definition such as the kinetic, the PS of the IS mass 
for phenomenological examinations. I will use the IS mass definition |HJ 
131 j for the discussions in the next section. The Greens functions also depend 
on the renormalization parameter v. For v of order |v|, the top quark 
velocity, the Greens function do not contain any large logarithmic terms 
and all logarithms are summed into the Wilson coefficients of the potentials 
and the currents by the renormalization group equations discussed in the 
previous section. Typically, one chooses v ~ 0.15—0.2, which corresponds to 
a momentum scale mtv ~ 25—35 GeV and an energy scale mtv 2 ~ 4— 7 GeV. 

Recently, in Ref. the order In a s corrections to the heavy quarko- 
nium partial width into a lepton pair were computed in the fixed-order 
expansion using an asymptotic expansion of QCD diagrams close to thresh- 
old |32| . (For a discussion of the behavior of the perturbative series in the 
fixed-order expansion I refer to Ref. 32 .) A summation of logarithms is 
not contained in that work, but the result can be used as a non-trivial 
cross check for the summations contained in the NNLL vector R-ratio R v 
of Eq. ()19|) . Expanding out the summations contained in the Wilson coeffi- 
cients for energies on the bound state poles one can determine the analogous 
order ajj lna s terms. The result agrees with the updated result of Ref. 32 . 
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0,0 t — , — i — , 1 1 1 1 1 1 1 1 1 1 1 1 — j 

346 347 348 349 350 351 352 353 354 




0,0 t — i 1 , 1 , 1 , 1 , 1 , 1 , 1 lJ 

346 347 348 349 350 351 352 353 354 

Fig. 3. Panel a) shows the results for Q 2 t R v with M 1S = 175 GeV and T t = 
1.43 GeV in fixed-order perturbation theory at LO (dotted lines), NLO (dashed 
lines) and NNLO (solid lines). Panel b) shows the results for Q\R? with the same 
parameters in renormalization group improved perturbation theory at LL (dotted 
lines), NLL (dashed lines) and NNLL (solid lines) order. For each order curves 
are plotted for v — 0.15, 0.20, and 0.3. The effects of initial state radiation, 
beamstrahlung and the beam energy spread at a e + e~ collider are not included. 



5. Discussion 

Let me now turn to how the NNLL non-mixing contributions affect 
the vector-current-induced top threshold cross section R v numerically. In 
Fig.Oh I have displayed Q%R V up to NNLL order. The curves show the LL 
(dotted blue lines), NLL (dashed green lines) and NNLL (solid red lines) 
cross section for v = 0.15,0.2 and 0.3. Figured shows the corresponding 
results in fixed-order perturbation theory. Compared to earlier analyses 
where the NNLL non-mixing contribution were not yet accounted for |121 
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I19| 120] . the NNLL cross section is shifted upwards by about +10%. On 
the other hand, the scale variation of the NNLL result is moderate and 
about ±3% for the variation of v I have used in this analysis. Compared to 
the fixed-order results with the same scales the improvement is substantial, 
particularly around the peak position and for smaller energies, but not as 
dramatic as concluded from our earlier analyses when the NNLL non-mixing 
contributions were not yet included. Although it appears premature to me 
to draw any definite conclusions from this result for the determination of the 
top Yukawa coupling, the strong coupling or the top quark width, because 
the mixing contributions are still not computed, it is prudent to say that 
presently the error estimate of ±3 % made in Refs. |121 1191 12(J| cannot be 
upheld and should be enlarged to b<Jttl tt ~ ±6% due to the relatively 
large shift between the NLL and the NNLL order results. As far as the 
determination of the top quark mass from a threshold scan is concerned, 
the new NNLL order results are not expected to affect the prospects of 
a determination with brat ~ 100 MeV since already for the fixed-order 
NNLO and the earlier NNLL order results the conclusions for the top mass 
determination from simulation studies were quite similar. [2J |3U] 



6. Conclusion 

In this talk I have discussed the ingredients needed to carry out a renor- 
malization group improved QCD computation of the total cross section 
cr(e + e~ — > tt) in the threshold regime. In renormalization group improved 
perturbation theory QCD logarithms of the top quark velocity are summed 
up to all order in a s according to the schematic expansion shown in Eq. ©. 
The summations can be obtained within an effective theory where all de- 
grees of freedom that can fluctuate close to their mass-shell are represented 
as fields and all off-shell fluctuations are integrated out. All logarithmic 
terms are associated to UV divergences in the effective theory and the sum- 
mation of logarithms is achieved by solving the renormalization group equa- 
tions of the couplings and coefficients of the effective theory. At present all 
ingredients for a full NNLL order prediction of the total threshold cross 
section are known except for the NNLL order result of c\, the coefficient 
of the dominant current that produces the top pair in e + e~ annihilation, 
which is only fully known at NLL order. I have presented new results for the 
NNLL non-mixing contributions of anomalous dimension of the coefficient 
c\. Numerically, the NNLL non-mixing contributions to c\ are of the same 
size as the NLL contributions and shift the total cross section by about 
+10%. Including the new results the present theoretical uncertainty for the 
normalization of the total cross section is about ±6%. The new result does 
not affect the prospects for a precise determination of the top quark mass 
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from a threshold scan based on earlier work, but it does affect the prospects 
for the determination of the top Yukawa coupling, the strong coupling and 
the top quark width. However, it is premature to draw definite conclusion 
as long as the NNLL mixing contributions have not yet been determined. 
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